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L. Grigoryan^ 

^Yerevan Physics Institute, Br.Alikhanian 2, 375036 Yerevan, Armenic^ 

Nuclear attenuation of the multi-hadron systems in the string model is considered. The im- 
proved two-scale model with set of parameters obtained recently for the single hadron attenuation is 
used for calculation of the multiplicity ratios of the one-, two- and three-hadron systems electropro- 
duced on nuclear and deuterium targets. The comparison of the features of the one-, two- and three- 
hadron systems is performed. The predictions of the model for multiplicity ratios of multi-hadron 
systems as functions of different convenient variables are presented. 
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I. INTRODUCTION 



In any hard process the initial interaction takes 
place between partons, which then turn into the final 
hadrons by means of hadronization process. However 
the hadronization process cannot be described in the 
framework of the existing theory of the strong interac- 
tions (perturbative QCD), because of major role of "soft" 
interactions. Therefore the experimental and theoretical 
(on the level of phenomenological models) studies of the 
all aspects of the transition from partons to hadrons are 
very important. The space-time evolution of the hadro- 
nization process, despite its importance, has been stud- 
ied relatively little. The study of the early stage of the 
hadronization process can shed additional light on the 
further development of the process. Semi-inclusive re- 
actions with nuclear targets give the possibility to study 
the development of the hadronization process on dis- 
tances of a few Fermi from the point of initial interac- 
tion. 

In particular, the nuclear attenuation (NA) of the 
high energy hadrons is the well known tool for inves- 
tigation a early stage of hadronization process ^ . There 
are many phenomenological models, which describe, 
rather qualitatively, existing experimental data for sin- 
gle hadron NA ||Il]-(n|]. Also some predictions for the 
attenuation of multi-hadron systems leptoproduced in 
nuclear matter in the framework of the string model 
were done |jl2]-(i3|]. It was argued that measurements 
of NA of a multi-hadron systems can remove some am- 
biguities in determination of the parameters describing 
strongly interacting systems at the early stage of parti- 
cle production: formation time of hadrons and cross- 
section for the intermediate state to interact inside the 
nucleus. Then, for the first time, data on the two-hadron 
system multiplicity ratio were obtained in electropro- 
duction llT3l . Experiment was performed in specific con- 



ditions. The multiplicity ratio of the charged hadrons 
was measured as a function of the fractional energy of 
the subleading hadron Z2, whereas over the fractional 
energy of the leading hadron zi the integration in the 
region 0.5 < zi < 1 — Z2 was performed. Later the data 
on the two-hadron system multiplicity ratio in neutrino- 
production were presented by another experiment [15]. 

The data on the two-hadron system multiplicity ra- 
tio [14] were described in the framework of some theo- 
retical models: the probabilistic coupled-channel trans- 
port model fl^, the so called energy loss model 1 17] and 
the string model [18]. In particular we showed [15, 18], 
that based on the two-scale model (TSM) [4] and im- 
proved two-scale model (ITSM) [10], it is possible to de- 
scribe these data quantitatively in the framework of the 
string model. We presented also predictions for the de- 
pendence of two-hadron system NA on the virtual pho- 
ton's energy in the same model. Possible mutual screen- 
ing of the hadrons occurred from one string and its ex- 
perimental verification have been discussed. 

In this work we continue the study the electropro- 
duction of multi-hadron systems in cold nuclear matter. 
This is the main goal of the present paper to consider 
the mutual screening of the prehadrons and hadrons in 
string (jet) ^. We compare one-, two- and three-hadron 
systems and show that mutual screening of prehadrons 
and hadrons plays essential role and can be measured 
experimentally. For instance such data can be obtained 
by HERMES Experiment, SKAT Experiment, and JLab 
after upgrade to the energy 12 GeV. We suppose that in- 
vestigation of the mutual screening of prehadrons and 
hadrons m cold nuclear matter can help to establish ini- 
tial conditions for the study of similar processes m hot 
nuclear matter arising in high energy hadron-nucleus 
and nucleus-nucleus interactions at RHIC and LHC. 

The paper is organized as follows. In Section II the 
theoretical framework is briefly described. Results and 
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^ NA means the difference of the ratio the multiplicities (per nucleon) 
on nucleus to that on deuterium from unity. 



^ The term "string" here means the object arising in result of DIS, 
which during its space-time evolution turns into states consisting 
of strings, prehadrons and hadrons. After all this object turns into 
the jet of hadrons. 
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FIG. 1: Leptoproduction of three-hadron system on nuclear target. 
Details see in text. 

discussion as well as necessary ingredients for calcula- 
tions are presented in Section III. Our conclusions are 
given in Section IV. 

II. THEORETICAL FRAMEWORK 

In the works (Tl-fTsI] the process of leptoproduc- 
tion of multi-hadron systems on a nucleus with atomic 
mass number A was considered theoretically for the first 
time. Although in these papers discussions were pre- 
sented for the general case of the n hadrons observed 
in the final state, and some formulae have been writ- 
ten for this general case, however the basic formulae 
and numerical calculations were performed for the case 
of two observed hadrons. In this work we make next 
step in this direction and consider three hadron sys- 
tems observed in the final state of the high energy semi- 
inclusive lepton-nucleus interaction. We do not give the 
equations for the case of one or two hadrons, which 
will also be used for calculations and discussions in the 
article because the conversion of these formulae to the 
case of one or two hadrons are very simple. The semi- 
inclusive reaction of the leptoproduction of three had- 
rons on nuclear target is: 

k + A^lf + hi + h2 + h'i+X, (1) 

where li{lf) is the initial (final) lepton, hi, /12 and are 
the observed hadrons. The hadrons hi, /12 and h^ carry 



fractions zi, zi and 23 of the total available energy (the 
energy conservation implies the condition: Z1 + Z2 + 23 < 
1). The multiplicity ratio for that process is defined as (it 
is assumed that averagings over transverse momenta of 
the final hadrons performed). 

= 2^^CTA(^^,(3^2l,^2,^3)/ 



Ada oiy, , ^'1,^2, Z'i) , (2) 

where duA and dan are the cross-sections for the reac- 
tion (1) on nuclear and deuterium targets, respectively, v 
denotes the energy of the virtual photon and — —q^, 
where is the square of the four-momentum of the vir- 
tual photon. One can imagine the reaction (1) as shown 
in Fig.l. The interaction of the lepton with the intranu- 
clear nucleon occurs at the point (6, x) from which the 
intermediate state q begins its propagation (6 and x are 
the impact parameter and the longitudinal coordinate of 
the DIS point). Initially the intermediate state q presents 
itself the object like a string with knocked-out quark on 
the fast and nucleon remnant on slow ends, which con- 
nected by means of string consisting of gluons. During 
further movement string breaks on the smaller pieces, 
and in result at the points (5, Xd), (h,Xc2) and \j3,Xr_z) 
the first constituents (valence quarks or antiquarks) of 
the hadrons hi, /12 and hj, are produced, and at the 
points (h,xi), {b,x2) and (6, ^3) the second constituents 
are produced and the yo-yo of the hadrons hi, /12 and ^13 
arise (the term "yo-yo" means, that the colorless system 
with valence contents and quantum numbers of the final 
hadron is formed, but without its "sea" partons). The 
points (6, Xci), {b, Xc2) and (6, Xcz) do not represented in 
figure, but they are used properly in the calculations. 
In Fig.l we for the sake of simplicity represent the case 
of three adjacent hadrons. In fact, all possibilities have 
been considered in the calculations both adjacent and 
not adjacent hadrons. 

We do not take into account the hadrons produced 
in result of decay the resonances. This factor could lead 
to an increase of nuclear attenuation, if taken into ac- 
count only one hadron from each resonance and a de- 
crease of nuclear attenuation, if taken into account that 
two or three hadrons can be produced from the same 
resonance. We think that the overall effect is small. 

In the string model there are simple connections be- 
tween above mentioned points xi—Xd ~ ZiL, X2—Xc2 = 
Z2L and xs — Xc3 — z^L, where L is the full hadronization 
length, L = v/ K, k is the string tension (k = IGeV / fm). 
The multiplicity ratio for the case of three hadrons ob- 
served in the final state can be presented in the form: 
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^ ^ J —OQ J X J xi J X2 

[D{zi,Z2,Z3,xi -x,X2 -x,X3 - x)Wo{hi,h2,h3;b,x,xi,X2,X3)+ 

D{zi,Z3, Z2,Xi - X,X2 - X,X3 - x)Wo{hi,h3, /12; 6, X, Xi,X2,X3) + 
D{Z2, Zi, Z3,Xi - X,X2 - X,X3 - x)Wo{h2,hi,h3] b, X, Xi,X2,X3) + 
D{Z2, Z3, Zi,Xl - X,X2 - X,X3 - x)Wo{h2,h3, hi] b, X, Xl,X2,X3) + 
D{Z3, Zi,Z2,Xi - X,X2 - X,X3 - x)Wo{h3, hi, h2;b, X, Xi, X2, X3) + 

D{z3, Z2, zi,xi - X, X2 -x,X3- x)Wo{h3, h2, hi;b, x, xi, X2, 2:3)], (3) 

where D{zi, Z2, Z3, li, I2, 13) (with li < I2 < I3) is the distribution of the constituent formation lengths li, I2 and I3 
of the hadrons and p{b, x) is the nuclear density function normalized to unity. Wq is the probability that neither the 
hadrons hi, /12, /13 nor intermediate states leading to their production (initial strings) interact inelastically in nuclear 
matter: 

Wo{hi,h2,h3\b,x,xi,X2,X3) ^ {I - Qi - {Hi +Q2+ H2 + Q3 + H3- 



Hl{Q2 +H2+Q3+H3- H2{Q3 + H3)) - H2{Q3 + H3))) 



(A-l) 



(4) 



where Qi, Q2 and Q3 are the probabilities for the initial strings of the corresponding hadrons to be absorbed in 
the nucleus within the intervals (x, xi), (xi, X2) and (x2, X3), respectively. Hi {i — 1, 2, 3) is the probability for the 
hi to interact inelastically in nuclear matter, starting from point Xi. The probabilities Qi, Q2, Q3, Hi,H2, H3 can be 
calculated using the general formulae: 



P(x, 



min 1 -^max ) 



app{b,x)dx , 



(5) 



where the subscript P denotes the particle (initial string or hadron), ap its inelastic cross section on nucleon target, 
and Xjnin arid x^ax are the end points of its path in the x direction, as it is shown in Fig.l. 

We use the scaling function of the standard Lund model for calculations. The simple form of this function 
f{z) ~ (1 + c)(l — zY, where c « 0.3 is the parameter which controls the steepness of the standard Lund model's 
fragmentation function, allows to sum the sequence of produced hadrons over all ranks and to obtain the analytic ex- 
pression for the any number of particles observed in final state. In the general case of the n hadrons the distribution 
D{zi, Zn] li, In) of the constituent formation lengths Zi,...,/„ is: 



D{zi---z„;li---L)^L''{l + cy 



ih--- InY 



iih + ziL) •••(/„ + ZnL)Y+- 



5{ln - (1 - Zn)L) + 



1 



5{ln - ln-1 - Zn-lL) + 



1+C 



+ Zn—lL 



S{l2 - h - ziL) + 



1 + C 

li + ziL 



(6) 



where Z„ < (1 — Zn)L, In-i < In — Zn-iL,..., < h < h — ziL. Equation (6) was obtained for the first time in Ref. 11311 . 
Unfortunately, corresponding equation (2.21) from Ref. 1, 13,1 contains some mistakes and uncertainties. 



III. RESULTS AND DISCUSSION 



example in the electroproduction process when in final 



Multi-hadron production depends on many vari- 
ables. This complicates the study of such systems. For 
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state are observed n hadrons, even after averaging 
over virtuality of photon and transverse momenta of 
hadrons the ratio of multiplicities depends on n + 1 
variables (the fractional energies of hadrons and energy 
of photon). Although we restrict ourself in this paper 
by three hadrons observed in final state, the simulta- 
neous consideration of four variables is very difficult 
especially in experimental study. We escape this dif- 
ficulty by inclusion of some additional averagings. 
We will consider following combinations of fixed and 
averaged variables: (i) the dependence on the fractional 
energy of one of the hadrons, the "trigger" hadron (z/^), 
whereas integrations are performed over the fractional 
energies of other hadrons and the energy of virtual 
photon u is kept at fixed value (in this paper it is fixed 
at value lOGeV); (ii) the dependence on the number 
of observed hadrons n, whereas the averagings are 
performed over ztr in some regions and u is kept fixed; 
(iii) the i/-dependence at fixed value of the "trigger" 
hadron fractional energy ztr = 0.3; (iv) the depen- 
dence on the fractional energy of multi-hadron system 
Z = J2^=i ^i' where Zi is the fractional energy of i-th 
hadron, n = 1, 2, 3 is the number of hadrons observed 
in the final state; (v) the dependence on n, where the 
fractional energies of the all observed hadrons are 
integrated in the region 0.1 < 2; < 0.33. 

At present it is assumed that hadrons produced from 
one string attenuate independently (full attenuation). 
This seems strange for the following reason. String has 
transverse dimensions comparable with the transverse 
size of the hadrons (at least no more). Therefore it is nat- 
ural to suppose, that hadrons produced from one string, 
may partially screen one another, what in result must to 
lead to the weakness of NA (partial attenuation). For 
the study of this effect, and for comparison with the 
basic supposition that hadrons attenuate independently 
(full attenuation), we consider partial attenuation in 
extreme case, when hadrons fully screen one another, 
and in result multi-hadron system attenuates as a single 
hadron. In accordance with above suppositions we con- 
sider four different cases for nuclear attenuation: (i) all 
parts of string and all produced hadrons are absorbed 
in nuclear medium independently (full attenuation). 
The full attenuation corresponds eq.(4). In all figures 
for notation of full attenuation we use solid lines; (ii) 
only initial string for the first produced hadron and first 
produced hadron itself attenuate (partial attenuation). 
Here the first produced hadron means the hadron first 
produced on time among observed ones. The partial 
attenuation corresponds eq.(4) with the corresponding 
replacements Q2 = H2 = Q3 = H3 = 0. In all figures for 
notation of partial attenuation we use dashed lines; (iii) 
it is supposed that n observed hadrons (n = 1,2,3) are 
adjacent ones and also that they produced on the fast 
end of the string. It is assumed additionally that this 




FIG. 2: Ratio R'^^'^^'^'^ for Krypton target at the energy of the vir- 
tual photon V = lOGe V as a function of the fractional energy of 
the "trigger" hadron ztr- Solid lines correspond to the case of full 
attenuation. From up to dawn the ratios Rm, Rm and R}^ are pre- 
sented, respectively. Dashed lines correspond to the case of partial 
attenuation. From up to down the ratios Ftf} and Rf} are presented, 
respectively. Dotted lines correspond to the case of full attenuation 
and adjacent hadrons produced on the fast end of the string. From 
up to down the ratios R]^, Rm and Rm are presented, respectively. 
Dot-dashed lines correspond to the case of partial attenuation and 
adjacent hadrons produced on the fast end of the string. From up to 
down the ratios Rm and R^ are presented, respectively. 



system suffers full attenuation in nuclear matter ^. The 
case of n adjacent hadrons produced on the fast end of 
the string corresponds eq.(6), where only (5-functions 
in square brackets are taken into account. In all figures 
for notation of n adjacent hadrons on the fast end of 
the string and full attenuation we use dotted lines; (iv) 
only n adjacent hadrons produced on the fast end of 
the string are taken into accoimt. It is supposed also 
that only initial string for the first produced hadron and 
first produced hadron itself attenuate. In all figures for 
notation of n adjacent hadrons on the fast end of the 
string and partial attenuation we use dot-dashed lines. 
Results of calculations with these conditions are 



' It is supposed that n hadrons observed in the final state are neigh- 
bors on the time of production, i.e. between them do not produced 
additional hadrons. The term "hadrons produced on the fast end 
of the string" means that it is a sequence of hadrons having lowest 
ranks in the string. 
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FIG. 3: Ratio Rlj for the Krypton target at the energy of the virtual 
photon V — IQGeV as a function ofn, where the fractional energy of 
the "trigger" hadron is averaged in the region: (a) 0.1 < ztr < 0.25; 
(b) 0.25 < Ztr < 0.4; (c) 0.4 < ztr < 0.55; and (d) 0.55 < 
Ztr < 0.7. The fractional energies of other hadrons are integrated 
over kinematically allowed regions. Notations are the same as Fig.2. 



shown in Figs.2-6. The nuclear density functions and 
set of parameters used in calculations were taken from 
our recent work I19ll. 



FIG. 4: Ratio Rfj'"^*"'^^ 



In Fig.2 the multiplicity ratios i?^^'^'''^'' for Krypton 
target at the energy of the virtual photon v — IQGeV 
as a function of the fractional energy of the "trigger" 
hadron ztr are presented. Solid lines correspond to 
the case of random selection of hadrons from the jet 
and full attenuation. From up to down the ratios R]^, 
and _R^J are presented, respectively. Dashed lines 
correspond to the case of random selection of hadrons 
from the jet and partial attenuation. From up to down 
the ratios and i?^'/ are presented, respectively. 
Dotted lines correspond to the case of adjacent hadrons 
produced on the fast end of the string and full attenua- 
tion. From up to down the ratios R]l}, Rf'} and i?^^ are 
presented, respectively. Dot-dashed lines correspond to 
the case of adjacent hadrons produced on the fast end 
of the string and partial attenuation. From up to down 
the ratios and Rf} are presented, respectively. 
From Fig.2 it is easy to see that high values of ztr 
{ztr > 0.7) are very convenient for the studying the 
mutual screening of the hadrons. For these values of 
Ztr positions of the particles in the string do not play a 
essential role. Mutual screening leads to the fact that the 
curves corresponding full and partial screening quite 
substantially different. The difference for the three- 



for Krypton target at fractional energy of the 
trigger" hadron ztr fixed at ztr ~ 0.3 as a function of the energy of 
the virtual photon v. Notations are the same as Fig.2. 
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FIG. 5: Ratio R^'^'^^'^'^ for Krypton target at energy of the virtual 
photon V = lOGeV as a function of Z = X]"-! ''^^^^^ 
fractional energy ofi-th hadron, n = 1,2, 3 is the number of hadrons 
observed in the final state. Notations are the same as Fig.2. 
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FIG. 6: Ratio at energy of virtual photon v = lOGe V as a 
function of n, where n = 1, 2, 3 is the number ofhadrons observed 
in the final state. The fractional energies of the all hadrons are inte- 
grated in the region 0.1 < z < 0.33. The results for different nuclei 
are presented: (a) Helium; (b) Neon; (c) Krypton and (d) Xenon. 
Notations are the same as Fig.2. 



particle case is more than for two-particle case. The 
positions of the particles in the string become significant 
in the case of small values of ztr {ztr < 0.3). Particles 
produced on the fast end of the string are attenuated 
less than others. The largest difference occurs in the 
case of single hadron, the smallest difference occurs in 
the case of three hadrons. 

In Fig.3 the ratio for the Krypton target at the 
energy of the virtual photon = lOGeV as a function 
of n are presented, where n = 1,2, 3 is the number 
of hadrons observed in the final state. The fractional 
energy of the "trigger" hadron is averaged in the 
region: (a) 0.1 < ztr < 0.25; (b) 0.25 < ztr < 0.4; 
(c) 0.4 < Ztr < 0.55; and (d) 0.55 < Ztr < 0.7. The 
fractional energies of other hadrons are integrated over 
kinematically allowed regions. Notations are the same 
as Fig.2. From Fig.3 we see that the study of small 
and medium ztr (panels a,b,c) in the case of n = 1 
can provide information about whether the observed 
hadron leading (i.e., that it was produced on the fast 
end of the string and contains knocked out quark) 
or not. In the case of n = 2 it is convenient to study 
positions of hadrons in string at small and medium ztr 
(panels a, b, c). In the case of n = 3 the study of medium 
and large ztr (panels c, d) can show the screening in the 
three-particle system takes place or not. 

hi Fig.4 the ratios i?^''^'''^'' for Krypton target at 



fractional energy of the "trigger" hadron ztr fixed at 
Ztr = 0.3 as a fimction of the energy of the virtual 
photon ;/ are presented. Notations are the same as Fig.2. 
Fig.4 shows that our questions: (i) attenuation is full or 
partial, (ii) observed particles are neighbors produced 
on the fast end of the string or not are, can be examined 
in the entire region considered energies. However, 
question (i) it is convenient to consider in the region 
relatively low energies {v ^ 5GeV), question (ii) in the 
region relatively high energies (v ~ 20GeV). 

In Fig.5 the ratios R^^''^^'^^ for Krypton target at 
energy of the virtual photon v = \QGeV as a function 
oi Z = Y^l^i Zi, where Zi is the fractional energy of i-th 
hadron are presented. Notations are the same as Fig.2. 
By this variable it is convenient to explore the question: 
are the produced particles neighbors produced on the 
fast end of the string or not? Behavior of the systems 
containing only the neighboring particles on the fast 
end of the string is qualitatively different from behavior 
of the systems that contain them among others. 

In Fig.6 the ratios RJ^ at energy of virtual photon 
V = lOGeV as a function of n are presented. The 
fractional energies of the all observed hadrons are 
integrated in the region 0.1 < ^ < 0.33. The results for 
different nuclei are presented: (a) Helium; (b) Neon; (c) 
Krypton and (d) Xenon. Notations are the same as Fig.2. 
It is easy to see that nuclear effects are amplified with 
the increasing of the atomic mass number A. The joint 
experimental study of the one-, two- and three-hadron 
systems can be very helpful. We want to note that in the 
energy region studied in this work (5-20GeV) the num- 
ber of hadrons in the current fragmentation region is 
limited and the probability to have among the observed 
hadrons two or even three adjacent ones is large. Also, 
it is likely that these hadrons were produced on the fast 
end of the string. 

Also we would like to briefly discuss why the cross 
section of the string-nucleon interaction may be equal 
to the cross section of hadron-nucleon interaction. Since 
the string is an object with small transverse dimensions 
the probability of mutual screening of the particles in 
the string is very large. 

There are other reasons why multi-hadron systems 
can attenuate as a single hadron. The two- or three- 
hadron systems will attenuate as a single hadron when 
final hadrons appear in result of decay of one resonance. 
For instance, combinations two or three pions can be 
obtained in result of decay of single vector meson 
produced in nucleus and decayed behind it. 



IV. CONCLUSIONS 

In this paper the problem of mutual screening of 
the prehadrons and hadrons in string in the framework 
of standard Lund model has been considered. We have 
shown that if the relevant data will be obtained, it will 
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assess the degree of mutual screening of the particles 
in the string. From our point of view, such informa- 
tion would be very useful for imderstanding the behav- 
ior of jets in high energy hadron-nucleus and nucleus- 
nucleus interactions. Unfortunately, many questions re- 
mained over the scope of this work: (i) How strongly 
the results depend on the chosen model? As mentioned 
above, a simple formula for scaling function in the stan- 
dard Lund model allows to obtain expressions for any 
number of hadrons in a compact form. Such compact 
expression can not be obtained in the case of more com- 
plex scaling functions (for instance, symmetric Lund 
model's scaling function), (ii) How the results change if 
we consider hadrons with specific charges and the dif- 
ferent cross-sections? We can choose a combination of 
particles that can not be neighbors in the string, or have 
a very different cross sections. Experimental study of 
such combinations can be very useful for the develop- 
ment of the model, (iii) How the results change if we 
consider that two or three hadrons could be produced 



as a result of decay of the single resonance? Above we 
tried qualitatively answer this question, but a quantita- 
tive study is needed, (iv) In this work the basic case was 
considered, when in the direction of virtual photon, in 
result of DIS, one string arises. In the work [20] we con- 
sidered the case, when in the direction of virtual pho- 
ton arise both one and two strings. Contribution of the 
events with two strings is relatively small in the case of 
single hadron. However it can essentially increase if in 
the final state are observed two or three hadrons. These 
questions will be discussed in further publications. 
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